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Abstract. Survival and percolation probabilities are most important quantities in 
the theory and in the application of growth models with spreading. We construct 
field theoretical expressions for these probabilities which are feasible for perturbation 
expansions. The outstanding role of the absorbing noise is stressed to obtain survival 
probabilities monotonic decreasing with time. We briefly consider some fundamental 
growth models equipped with absorbing noise which are representations of known 
universality classes of spreading phenomena. The critical scaling properties of their 
survival and percolation probabilities are stated. In an appendix we consider shortly 
the renormalized field theory of compact directed percolation. 
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1. Introduction 

The investigation of the formation and properties of random structures has been an 
exciting topic in statistical physics for many years (see e.g. [I1I21)- Of particular interest 
are random structures which are formed by local rules. These processes can often 
be expressed in the language of population growth. Stochastic processes of this type 
describe the essential features of a vast number of growth phenomena of populations 
near their extinction threshold, and are relevant to a wide range of models in physics, 
chemistry, biology, and sociology IH El El IZ| ■ The transition between survival and 
extinction of such a growing population constitutes a nonequilibrium continuous phase 
transition phenomenon and is characterized by universal scaling laws. The extinct states 
without active individuals (interacting particles) are stochastically absorbing, i.e. in a 
finite system with a finite number of degrees of freedom an extinct inactive state is 
reached with probability one at a finite time. Beyond the critical transition point, 
where the replication rate is greater than a critical value, the systems are in the active 
state. The corresponding mean extinction time is practically infinite in systems with a 
great number of degrees of freedom. To find active steady states it is therefore allowed 
to take the infinite volume limit, the so-called thermodynamic limit, before the infinite 
time limit is done. Then, the probability Poo that an active state, created by a localized 
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seed is still alive after an infinite time is greater than zero beyond the critical point. 
Hence, this percolation probability serves besides the density of the active particles 
in the steady state, pst, as an order parameter for growth processes. On the other 
hand, in the absorbing phase, the phase below the critical replication rate, each active 
configuration dies after a finite time. Thus Poo is zero here and the steady states are 
inactive. However, an interesting observable in both phases is the survival probability 
P{t). It is the probability that a process started from a localized seed, e.g. an active 
particle at the origin of space and time, is still active at time t. Of course, we have 
generally P{t) P^a for t — > oo. 

Near a critical point one is typically interested in the asymptotic universal long- 
wavelength, long-time behavior of the observables. Besides the existing vast number of 
simulation methods and results (see e.g. 0), a variety of analytical tools were developed 
for growth processes in the past decades. The method which introduces the continuum 
limit from the beginning and exploits the renormalization group to obtain universal 
properties is the renormalized stochastic field theory |Hl El HHl HI]- At the heart of 
this method is the calculation of time dependent mean values, correlation and response 
functions of the observables by using path integrals. Basic ingredient is a generalized 
probability measure exp(— ^J), where the dynamic functional J (nowadays also called a 
response functional) is a functional of time and space dependent fields. The restriction 
of J to the only relevant terms in an expansion of the slowly varying fields defines 
a renormalizable field theory. This theory represents then a full universality class of 
systems with common asymptotic properties. 

Although the field theory of growth processes are based on a probability measure, 
the direct calculation of the survival probability P{t) and the percolation probability Poo 
for growth processes by field theoretic means was not well developed in the past. There 
exists an attempt of Munoz, Grinstein, and Tu [12] (henceforth cited as MGT) to devise a 
practical method to determine correctly the scaling properties of the survival probability. 
However, these authors erroneously assert: "In field theories with continuous variables 
. . . the absorbing state is a set of measure zero in phase space and so can never actually 
be reached in finite time. Thus P{t) is strictly equal to unity for all t, so the concept of 
the survival probability has no utility." , an argument that resembles Zenon's paradox on 
Achilles and the tortoise fHI- They introduce therefore another quantity as a substitute 
for P{t) with the correct scaling properties. But their argumentation seems to rest solely 
on a naive interpretation of the Langevin equation as an ordinary differential equation for 
a continuously differentiable function, and of a stationary solution of the corresponding 
Fokker-Planck equation for the probability density of the stochastic process in zero 
spatial dimension. However, this stationary solution has a non-integrable singularity at 
= 0, and is therefore not normalizable. Hence, it cannot represent any features of the 
correct probability density near the absorbing state = 0. MGT, without considering 
any boundary condition of the time dependent solution of the Fokker-Planck equation, 
conclude: "An arbitrary initial probability distribution therefore evolves in time toward 
a distribution weighted at values lying progressively closer to zero. Note, however, that 
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cannot actually achieve the value in finite time." However, to reveal such properties 
of the time dependent probability distribution one has to ascertain analytically the 
behavior of the solutions of the Fokker-Planck equation near and at the absorptive state. 
But this was not done by MGT. Hence, their arguments are far of being convincing. 
We will show that the validity of their conclusion depends heavily on the form of the 
underlying noise. For a noise called absorptive in the following, the absorbing state has 
a 5-probability measure in the continuous phase space. This measure accumulates more 
and more mass during the course of time. Consequently, the survival probability P{t) 
is a monotonically decreasing function also in a continuous phase space. In contrast, for 
the so-called multiplicative noise jj^ the assertions MGT are correct, and P{t) = 1 for 
all finite times. 

In this paper we develop expressions for survival and percolation probabilities of 
several react ion- diffusion processes with absorbing states in non-zero spatial dimensions 
from which the calculation of these quantities by field-theoretical methods can be 
performed. We show the feasibility of our expressions for the calculation of nontrivial 
probabilities in a mean field (saddle point) approximation that becomes correct for 
spatial extended systems above their higher critical dimension. 

The outline of our paper is as follows. In section II we present a simple toy model 
that shows the fundamental difference between absorptive and multiplicative noise. It 
provides a nontrivial survival probability in the absorbing case. Section HI reviews well 
known field theoretic results for directed percolation (DP) and its properties under the 
duality transformation. In section IV we derive the field theoretic expressions of the 
survival and the percolation probability for DP. Asymptotic scaling forms are found that 
are related to the mean particle density of the dual process. In section V we calculate 
the probabilities in a saddle point approximation. In section VI we study a variant of DP 
equipped with multiplicative noise. We show that in this case the survival probability 
is equal to one for all finite times. We compare in this section our result for the survival 
probability with the expression introduced by MGT. In Sec. VII we apply our results to 
the field theories of various fundamental growth processes. A short epilogue is given in 
section VIII. An appendix contains a brief presentation of the renormalized field theory 
of compact directed percolation. 

2. A toy model 

For the ease of the reader we re-consider in this section a well known example of a 
continuous growth process for which survival probabilities can be easily discussed in a 
concise way. Most of this material can be found in textbooks [13 El CZl EHj , but it does 
not seem much appreciated by those working in this research field. The stochastic version 
of the logistic or Malthus-Pearl-Verhulst equation, x = —r{x)x + C with r{x) = (r + gx), 
describes the growth of a population size x{t) > 0. Equipped with absorptive noise (, it 
provides a zero- dimensional toy-model of directed percolation. The model is also useful 
for the mean field description of the evolution of a homogeneously distributed population 
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in a spatial extended system with volume V. In this case one defines the particle density 
by n{t) = x{t)/V. The coupling constant scales with the volume as (7 ~ 1/^- The 
Gaussian random noise ({t) is characterized by correlations C{t)C{t') = 2x{t)"6{t — t'). 
For physical growth processes these correlations as functions of x, t, and t' result from 
the elimination of many fast microscopic degrees of freedom in local equilibrium with 
the population size x. Therefore, they should be analytic functions in x and zero if 
X = 0. Hence, in the limit where microscopic correlation times go to zero on the slow 
macroscopic time-scale, we typically find for population processes a = 1 (absorptive 
noise) or a = 2 (multiplicative noise) in a small x expansion neglecting higher orders. 
The nowadays popular terminology "multiplicative noise" stems from the fact that in 
the before-mentioned time limit, the noise may be formally written as a product of an 
analytic function of x and of white noise. In any event, it is essential that the noise is 
interpreted in the Ito-sense I18j. A very qualitative glance on the logistic Langevin- 
equation shows that the random impulses resulting from the noise scale like Thus, 
for small x, these impulses can easily surmount the distance x to the absorbing state 
X = if a < 2. We expect that this is the condition for which a non-zero extinction 
probability Po{t) = 1 — P{t) for finite times can be defined. 

Nevertheless, a more satisfactory tool for the analytical description of this stochastic 
process is given by the corresponding Fokker-Planck equation for the probability density 
p{x, t) to find a population with x > at time t: 



The Fokker-Planck equation is also known as the forward Kolmogorov equation for 
the fundamental probability density p{x,t) = w{x,t\x' ,t') with t > t' and the initial 
condition w{x, t'\x', t') = 6{x—x') with < x' < 00. As a second order partial differential 
equation it needs in general two boundary conditions at x = and x = 00. Starting in 
the open interval (0, 00), it may be that the stochastic process x{t) does never takes on 
the value x = or x = 00 in finite times. In this case the corresponding boundary is 
called natural, and no boundary condition at this value has to be imposed. In contrast, 
if the process x(t) has positive probability of taking the value x = or x = 00, the 
corresponding boundary is called accessible. For accessible boundaries there are two 
possible cases: 

(i) Exit: the drift towards the boundary is such that the boundary automatically acts 
as an absorbing barrier, and no boundary condition can be imposed. 

(ii) Regular: the process behaves like a classical diffusion process near the boundary, 
and various boundary conditions can be imposed. 

Hence, for the existence of an unbiased nontrivial extinction probability, the absorbing 
state at X = should be an exit boundary. 

According to the scheme of Feller ^Sj, the classification of boundaries depends on 



dp{x,t) d[r{x)xp{x,t)] ^ d'^[x' 
dt dx 



dx"^ 



'■p{x,t)] 



(2.1) 
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the Lebesgue integrability of the function 




(2.2) 



where xq G (0, oo). The following functions are defined besides F{x): 



(2.3a) 



H 



(x) = F(x)-i / dyG{y) 




XQ 



{2.3b) 



Note that the stationary solution of the Fokker-Planck equation ()2.H1 . G{x), is a 
stationary probability distribution only if it is normalizable. We denote with f{x) G C{0) 
or f{x) E C{oo) a function that is Lebesgue integrable on the interval (0, Xq) or (xq, oo), 
respectively. Then Feller's classification criteria are: 

(i) The boundary b is regular if F{x)^^ G C{b), G{x) G C{b), and H{x) G C{b); 

(ii) The boundary b is exit if F{x)-^ G C{b), G{x) ^ C{b) and H{x) G £(6); 

(iii) The boundary b is natural in the other cases. 

Applying these criteria one finds easily that for a < 3 the boundary at infinity is natural, 
whereas the boundary at zero is regular for a < 1, exit for 1 < a < 2, and natural for 
2 < a (attractive for r > —1, repulsive for r < —1). It is now easily seen that for 
1 < a < 2 the Fokker-Planck equation p.lj) is transformed back to the case a = 1 by 
the setting y = and the replacement r{x) —>■ r{y) = r{x) + (a — l)/(2 — a). Hence, 
in the following we concentrate on the absorptive noise a = 1 and to the multiplicative 
noise a = 2. In both cases no additional boundary conditions are needed for the solution 
of the initial problem. 

On the other hand we find from above that there is a fundamental difference 
between these two cases: for absorptive noise there is a non-zero probability for the 
stochastic process to reach the absorbing state at a; = in finite times, whereas 
for multiplicative noise the probability is zero that the process reaches x = in 
finite times. This different behavior is reflected in the evolution of the probability 
distribution in time: for absorptive noise the probability density shows a 5-peak directly 
at the absorbing boundary which accumulates more and more mass in the course 
of time. For multiplicative noise in the case r > — 1 the probability distribution 
p{x,t) represents a peak, which moves to the left, becomes smaller and smaller in 
extension, and accumulates all probability near the absorbing state x = without 
reaching this state itself {p{0,t) = for t < oo). Only in the infinite time limit we get 
limt-,oo p{x , t) = 6{x) and the stochastic process is ultimately absorbed. For r < — 1 
the probabihty distribution with p{0,t) = for t < oo converges to the stationary 
distribution limt^oo p{x,t) oc G{x), equation ()2.3aj) . However, note that the limits 
limt_^oo and lim^-^o cannot be interchanged |2Uj . 
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The fundamentally different behavior of the distributions in time can be made fully 
explicit in the linear case r{x) = r [THl ITHl 1^. If the noise is absorptive, a = 1, the 
solution of the Fokker-Planck equation ^2.\\ is found easily by a Laplace transform to 

t\xQ) = 5{x)Po{t, Xo) + 9{x)p{x, t\xo) , (2.4a) 
pix,t\xo) = /i(t)Vi(t)7a;/i(2/i(t)v/i(t)^)exp[-/i(t)(x(t) +x)] , (2.45) 
Po{t,Xo) = exp[-h{t)x{t)] . (2.4c) 

Ii is the modified Bessel function of first order, and we have defined 

x{t) = Xo exp(-rt) , {2Ad) 

hit) = . (2.4e) 

1 — exp(— rt) 

Note that the extinction probability (j2.4cp increases in the course of time monotonically 
from to 1 (or exp(— |r| xq) if r < 0, whereby Poo(a^o) = linit^oo(l — -^0(^52:0)) = 
1 — exp(— |r| a;o) is a kind of percolation probability). Hence, it is explicitly shown that 
the assertion "the stochastic process cannot actually achieve the value in finite time" 
of MGT is wrong in the case of absorptive noise. 

The behavior of the probability distribution is qualitatively different in the case of 
multiplicative noise, a = 2. Then, the Fokker-Planck equation ()2.H) is solved by the 
substitution x = exp y that transforms the stochastic process to simple diffusion with 
drift. One obtains 

p{x,t\xo) = {V^x)'^exp - [lna;/a;o + {T + l)tf/At^ . (2.5) 

Indeed, here p{0,t\xo) = and Po(t,xo) = for all finite times, and the assertion of 
MGT is correct. 

The fundamental probability density w{x,t\x' ,t') must also be a solution of the 
backward Kolmogorov equation 

dw{x,t\x',t') ,dw{x,t\x',t') ,^d^w{x,t\x',t') 

m' = ^^"^" — d^' " — 87^ — • ^'-'^ 

With the Ansatz w{x,t\x' ,t') = 6{x)Po(t — t',x') + 9{x)p{x,t\x' ,t') one formally 
extends the range of the variable x in ()2.1|) to the open interval (—00,00). Then 
dxp{x, t\x', t') = Po{t — t', x') + dxp{x, t\x', t') = 1 where p{x, t\x', t') fulfills the 
original Kolmogorov equations. For Po{t,x') one gets the ordinary differential equation 

d r°° 

Po{t-t',x') = -— dxp{x,t\x',t') = alim{x''-^p{x,t\x',t')) . (2.7) 
dt Jo x^o 

Po is zero for multiplicative noise but equal to p{0, t\x', t') for absorptive noise. From the 
backward equation ()2.fij) on derives the differential equation of the survival probability 
P(t, x') = 1 — Po(t\x', t' = 0) starting from a state with population size x' > at t' = 0. 
Using the time translation invariance of the process, the time evolution follows from the 
backward Kolmogorov equation ()2.6p as 

dP(t,x) , , dPit.x) ^d^P(t,x) 

^ = (2.8) 
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with P(0,x) = 1 for X > 0. Using this equation the mean time T{x) = dt P{x,t) to 
reach the absorbing boundary at x = from a state a; > is easily calculated 



n^)= dy ^^exp(-/ dzz'-'^riz)). (2.9) 



ds 

(sTy)" 

With r{x) = T + gx one obtains in the case of multiplicative noise 

PX POO 

^^''^ = I "^yyj^ (FTW^r^M-gs). (2.10) 

This expression diverges at the lower integration boundary, y = 0, irrespectively of the 
sign of r for all x > 0. Hence, as we already know, the stochastic process cannot reach 
the absorbing state in finite time. 

In the case of absorptive noise we obtain 
r /"°° ds 

T(x)=/ dy ^exp[-gsy2-{T + gy)s]. (2.11) 
Jo Jo s + y 

This expression is finite for all x > and shows again that the extinction probability 
to reach the absorbing state is nonzero even for finite times. From the mean field point 
of view, starting with a homogeneous state in the volume V , we have x ~ V and 
g ~ '^/V. Consequently, for r > the mean extinction time scales with the volume as 
T(x) ~ |lnx|, whereas for r < we have an exponential increase InT(x) ~ 'r'^/g- 
This profound difference between the absorbing and the active phase of population 
growth is also expected in extended but finite systems. 

Our discussion of known results from the analysis of continuous stochastic processes 
has shown that the behavior of the survival probability in time depends crucially on the 
type of the noise. It is shown that there exists a non-trivial survival probability in a 
theory with absorptive noise also in a phase space of a continuous variable. Thus, the 
toy model does not provide any argument against the existence of a non-trivial survival 
probability in the case of field theories of growth processes in spatial dimensions greater 
than zero. Indeed, in the next chapters we directly construct calculable expressions of 
survival and percolation probabilities for such field theories in terms of path integrals. 



3. Directed percolation 

In order to be consistent with later parts of the article, we will concisely present 
well known material on critical stochastic processes with absorbing states in this 
section. We will concentrate on directed percolation (DP) as the typical growth model 
ini 122 123 121 123 123 121 ■ For an comprehensive overview see E]. 

From the onset we use a mesoscopic picture in which all microscopic length- and 
time-scales are considered to be very short. We are interested in the asymptotic 
properties of the growth process near the critical point between the absorbing phase 
and the steady active phase. Thus, we take a continuum approach with the density 
of the active particles, n (x, t) , as the stochastic variables. We remove all irrelevant 
contributions to the equations of motion. The growth process might be represented by 
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a Langevin equation (in the Ito sense), which is constructed in accordance with the 
principle of the existence of an absorbing state as [23 |2Z1 12H] 

dtn = W^n- R[n]n + Q + C, (3.1) 

where C(x, t) denotes the noise which has to vanish if n(x, t) = 0. The reaction rate 
R [n] models birth and death of particles, as well as the saturation of the population. In 
addition we introduce an external particle source Q(x, t). The Langevin source Ci^^t) 
can be assumed to be Gaussian with absorbing correlation 



C(x, t)C(x', t') = Xg"n{^, t)5(x - x')5(t - t'). (3.2) 

Diffusive noise is irrelevant in comparison to the reaction noise, equation ()3.2j) . in the 
case of DP, but in general indispensable for processes with multiplicative noise [23 123- 
The extinction rate is given by 

R[n] = x(t + ^n) . (3.3) 



2 

Under the influence of fluctuations, the critical point is found at a value t = Tc < 0. We 
implicitly renormalize this parameter by r ^ r — Tc so that the critical point is always 
found at r = 0. 

The field (^(x, t) > of the Langevin equation (|3.ip describes the external source 
distribution of particles. As a special case the creation of one particle as a seed at the 
origin x = and t = is represented by a source term 

Q(x,t) = 5(x)5(t). (3.4) 

On the other hand, the choice Q = Xh = const describes the creation of particles at a 
constant rate Q uniform in space and time. In the following we are interested in the 
particular source model 

g(x, t) = Xh + po6{t) = X{h + po6{Xt)) , (3.5) 

where po constitutes a homogeneous particle density as an initial condition. 

Dynamic response functionals ^7 |Hl IHl QHI based on Langevin equations of the form 
()3.1|1 with noise correlation ()3.2j) have the general form [23 12S1 12H] 

J[fi,n]= y"dtd'^x|An(^A-i^ + (r- V^) + V[n,n])n|. (3.6) 

Here, n(x, t) denotes the response field conjugated to the particle density n{-x,t). For 
DP we have 

V[n,n] = ^{g'n-g"h). (3.7) 

The introduction of the source Q leads to the shift 

J[n, n] jQ[h, n] = J[n, n] - (Q, h) (3.8) 

where {Q, h) denotes the integral f dt d'^xQ{'x, t)n(x, t). The responses are defined with 
respect to the particle source Q(x, t) > 0. They are represented by correlation functions 
with the response field n(x, t). Such response and correlation functions (called Green's 
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functions in the following) can now be calculated by path integrals with the probability 
weight exp(— JT"). We use the notation 

{A[n,n]) = J V{n,n)A[n,h]exp{-J[h,n]), (3.9) 

where .4[n, n] is any functional of the fields n,n, and V{n,n) denotes the integral 
measure which is properly normalized (1) = 1. In a suitable discretization it reads 
V{h,n) = Y[y:ti^n{'x,t)dn{yi,t)/27ii). For practical calculations, the observables ^[n, n] 
should be either polynomial or exponential. Then a diagrammatic perturbation theory 
can be developed. The integration of the response fields is done along the imaginary 
axis, but it can be deformed to general complex values in finite regions. The integration 
of the density field can be formally extended over the full real axis. Initial and final 
conditions are supplied by 

n(x, — oo) = n(x, oo) = (3.10) 
for all X. For Q = we have therefore 

(n(x,t)) = (n(x,t)) = (3.11) 
for all x and t. 

Now we define the duality transformation (known as rapidity reversal in Reggeon 
field theory and corresponding to duality in the mathematical theory of interacting 
particle systems [Hj) as 

a'-^hi^x, t) =: 5(x, t) < > — s(x, —t) := — Q;n(x, —t) . (3-12) 

The free, bilinear, part of the dynamic functional ()3.(j|l is invariant under this 
transformation for each finite a. Note that has the dimension of a spatial volume. 
Choosing a = ^/g'/g" the interaction part of DP, V[fi, ra] (j3.7|) . is (up to the time 
inversion) also invariant under duality. Of course, a is a redundant variable from 
the renormalization group point of view. Such variables should be eliminated before 
any renormalization group consideration is applied. Choosing therefore s and s as the 
fundamental fields, the invariant coupling is g = \/g'g"i and the fields s and s have the 
same scaling dimensions in DP. For general growth models, the duality transformation 
defines a dual stochastic process which can be different from the original one. Such a 
growth model is therefore not self-dual like DP. The quantities belonging to the dual 
process are designed in the following by a hat (e.g. the mean particle density of the dual 
process is denoted by p). 

As a last point we note the asymptotic scaling properties of the Green's functions, 
the cumulants of the variables n and n. One obtains 

N N 

(nn(x.,t.)n-(x„t.))^^"'^^=irr^^^^ 

i=i j=i 

X i^;v.^({|rr X, t}, V \rf^''^'''' , Pol Vt^') (3.13) 

with universal functions -F/vtv- Particular scaling laws follow from this asymptotic result 
of the renormalization group equation. It shows that in general the scaling of growth 



Survival and percolation probabilities 



10 



processes is defined by four independent scaling exponents 2,^,13, j3 PHI- However, for 
theories wfiich are self-dual like DP we have in addition the relation (3 = (3. 

4. Percolation and survival probabilities for directed percolation 

The expressions for the percolation and the survival probability for DP which are 
constructed in this section, can also be used, mutatis mutandis^ for other growth 
processes, briefly considered in a later section. 

Percolative spreading processes with absorbing states like DP provide two 
independent fundamental order parameters, the density of the active particles in the 
steady state pst, and the percolation probability Poo- Near the transition point they 
scale as pst ~ and Poo ~ kl^ where, in general, the two exponents j3 and j3' are 
distinct from each other jHOj- It is well known that this asymptotic scaling behavior 
follows from the relations pst ~ {n) and Poo ~ (^), where, in particular, one has to 
specify the meaning of the expectation value for the latter expectation value. Due to 
the duality transformation ()H.12|1 . Poo is intimately related to the particle density of the 
dual process: Poo ~ Pst, so that (3' = (3 = (3. Hence, we get for self dual processes (like 
DP) the well known relation (3' = (3. In order to calculate scaling functions and universal 
amplitude ratios it is, however, not sufficient to consider only the asymptotic scaling 
behavior. We have to derive the exact relation between the two order parameters. 

The discussion of the cluster probabilities is complicated by the fact that one has 
to deal carefully with the probabilities coming from the non-scaling small clusters, 
say, with particle numbers smaller than Nq. For these small clusters the asymptotic 
continuum field theory cannot be used. However, the small clusters do not induce 
critical singularities due to their finiteness. We assume as usual that they provide only 
analytical corrections P]. 

Let Pn{t) be the probability for a directed cluster in [d + 1) -dimensional space- 
time consisting of particles created from one particle, the seed, at the origin of the 
space-time lattice, which is described by the source fj3.4|) . For N > \ one can use 

the continuum approximation 



The average in equation ()4.1|) is performed with the response functional, equation ()3.6|) . 
One expects the asymptotic scaling behavior P/v(t) = /{rNy) pp. In many cases, 
e.g. absorptive noise as in DP near the critical point, and ^ 1, it is possible to 
expand the exponential exp(n(0, 0)) in equation ()4.H) in a power series to get asymptotic 
results. Equation (jH.lOj) leads to {F[n]) = F[0] for any functional F of n. Hence, the 



PNir) ^ Pjv(r) = {5iN - M) exp(n(0, 0))) 



(4.1) 



where we have defined 




(4.2a) 



(4.2&) 
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first expansion term of tlie exponential, tlie 1, inserted into (|4.H) in place of expn(0, 0) 
leads exactly to S{N) = if ^ 1. The higher powers of n(0, 0) contribute non- 
leading corrections to the asymptotic scaling behavior. Hence, only the term n(0, 0) of 
the expansion must be retained. 

We define a generating function 

oo 

F(r,fc) = ^p^(r)e-'=^, k>0. (4.3) 

N=0 

This function has the property 

T^/ r ^\ 1 1 for r > , ,, 

^'^•*-'+°)=|l-P. for r^o: 

From the continuous counterpart of ()4.3|1 we obtain, neglecting possible analytic 
corrections that arise from the small clusters 

POO 

F{t, k) - F{\r\ ,k)= dN (P^(r) - PAr(|r|))e-'=^ 

^0 

= (exp(n(O,0) - A;7V))(r) - (exp(n(O,0) - A;7V))(|r|) . (4.5) 

Here we rely on the property that divergencies like Pjv ~ N^^ for small are 
independent of r and cancelled in the integrand of equation ()4.5|1 p. Using 
now the fact that in the subcritical absorbing phase with unbroken symmetries 
limfc^+o(e"^-^expfi(O,0))(r > 0) = (exp ri,(0, 0))(r > 0) = 1, we arrive at 

Poo = 1 - lim (exp(n(0, 0) - k^^)) . (4.6) 

Provided that the asymptotic expansion of the exponential is allowed we get 

Poo = - lim (e-'^^n(O,0)) oc Wf (4.7) 

by applying equation ()3.13|1 . With the help of the duality transformation ()3.12|1 in 
conjunction with the time translation invariance, we obtain finally 

Poo = a' lim (n(0, oo)e^'^^)('^'^"') = a^pst oc Irf (4.8) 

where J\f = f dt f d'^a;n(x, t) and Xh = k/a^ constitutes a constant particle source that 
is taken to zero at the end of the calculation. Note that the particle density of the 
dual process pst is independent of any initial state unequal to the vacuum. In self-dual 
theories like DP, (3 = (3 and the parameter a can be chosen that = Pst, as shown 
above. Then, we see that the redundant variable a enters the relation between the two 
order parameters: Pqo = d'^Pst 

Now we turn our attention to time dependent observables. As above, there are two 
fundamental complementary quantities, namely the particle density p{t) starting from a 
fully occupied initial state at time t = 0, and the survival probabihty P{t) that a cluster 
grown from a single seed is still active at time t. Both are monotonically decreasing 
functions of time. The exact relation between the two is found with the help of the 
following consideration. 
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Let PN{t) denote the probability that a cluster has exactly N active particles at 
time t. We define the function 

oo 

P{k,t) = Y,i^-e-''')p^{t). (4.9) 

N=0 

Of course, we have limk^oo P{k,t) = P(t). We distinguish again between small and 
large clusters. The former lead to uncritical corrections. For the large clusters we use 
the approximation 

PN{t)^PN{t) = {5{N-Ar{t))expn{0,0)) (4.10) 

where A/'(t) is defined in equation ()4.26|) . Thus, we get from the continuous counterpart 
of equation ()4.9j) 

poo 

Pik,t)= / diV(l -e-'=^)P^(t) = (expn(O,0)) - (exp(n(O,0) - A;7V'(t))) . (4.11) 
Jo 

The first term of the last equation, (expn(O,0)), is equal to one. Hence, we obtain the 
survival probability 

P{t) = 1 - lim (exp(n(0, 0) - kj^{t))) . (4.12) 

After expanding the exponential to get asymptotic results, and utilizing the time 
translation invariance of the process we arrive at 

P{t) = - lim (n(0, -t)e-'=^(°)) oc H/"" . (4.13) 

k—*oo 

Using again the duality transformation ()3.12|1 . we find 

P{t) = a" lim (n(0,t)e^o^(°))(^""i) _ ^2^^^^ ^ ^-p/vz ^ 

where M{Q) = J d'^xn(x, 0), and po = k/a"^ constitutes an initial constant particle 
density. Hence, p{t) is the particle density of the dual process at time t starting from a 
infinite homogeneous initial density po, corresponding to a fully occupied initial state. 
For self-dual processes like DP we have in fact p{t) = p{t) and j3/vz = f3/vz. Relation 
()4.14|1 is fully analogous to the steady state relation ()4.8|1 . but note the different limits 
/i — ^ and Po oo. In any case it follows from limt_^oo p(t) = Pst that 

lim Pit) = Poo . (4.14) 

t— »oo 

The expressions fOl 1^ KTB, KT^ and the relations (jlSl I^T^ are the main 
general results of this paper. They hold near criticality and for not too small times t, 
because we have neglected several analytic corrections and used, for discrete processes, 
the continuum approximation. Also we have neglected in ()4.7[ I4.13|) higher powers of n. 
Being interested in non-asymptotic results one has to use the full expression [expn — 1) 
instead of the single term n (see ()4.6|) and ()4.12|) ). 
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5. Probabilities in a mean field approximation 

To demonstrate the applicability of our results ()4. 7^14. 13]) . we calculate the corresponding 
path-integrals in a mean-field approximation. Then we are able to write down explicit 
simple expressions for all quantities and relationships. At first we determine the 
asymptotic expression of the cluster probability ()4.1|1 . Instead of directly evaluating 
this quantity we consider the Laplace transform 

POO 

M{k,T) = / dNe-^^PN{T) = (exp(n(O,0) - kM)) 
Jo 

V{h, n) exp{-J[h, n] - {k, n) + n{0, 0)) . (5.1) 

The saddle point of the exponential weight follows from the solutions of the variational 
equations 

5(J[h,n\ + (k,n)) dn „2\~ a / ~ ^q" ~2 , « ^ \ 

dn{x,t) ot 2 

6(J'\n,n] + (k,n)) dn Xq' o , „ _ , 

a(x.«) = a + - ^ + i" - ^» = ° • <^-2*' 

Replacing k by k6{T — t) where the limit T — oo is performed at the end, the solutions 
of the saddle point equations, subject to the constraints ()3.10|) . are 
n(xt)| -mT-t){eMMT-t))-l)k 

n(x, t)U,,, - ^ exp(A«;(T - t)) + {w - r) ' ^^'^"^ 
n(x,t)|,.p. = 0, (5.35) 

where w = a/t^ + 2g"k/X. Consequently we have (J7'[n, ra] + {k,n))\s.p. = 0, and we 
obtain in the limit T —>■ oo 

M{k,T) = expn(O,0)|,.p. = exp (^(r - + 2g"k/X)/g"^ . (5.4) 

Using equation ()4.6p . the limit k +0 of this expression yields the percolation 
probability 

P^{t) = 1 - hm M{k, r) = 1 - exp((r - |r|)//) ^ ^ ^(-r) . (5.5) 

The last expression is in accordance with equation (j4.8|) . since Pst(r) = 2 |r| 9{—T)/g' 
is the stationary solution of the mean field equation. 

From equation ()5.Hl we find the cluster probability P/v by an inverse Laplace 
transform of M{k, r) 

PN{r) = — / dke''^M{k,T) = -j—r exp ( - —{Xt'N+{XN)-'] 
Ar-3/2 

^ -=== exp{r/g" - Xr^N/2g") (5.6) 

y/2T[Xg" 

The last row presents indeed the well known asymptotic probability distribution for the 
generation of a directed iV-cluster on the Bethe-lattice p. Pnij) shows the symmetry 
property 

^A^(-kl) = (l-i^oo(-|r|))P^(|r|) (5.7) 
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between the probability distributions above and below the percolation point with 

rex 

/o 



dN P^i\r\) = 1. 



Now we turn our attention to the survival probability (j4.12j) . We have to solve the 
saddle-point equations ()5.2aj) where k is replaced by k6{t) or, in other words, with the 
final-condition ?7,(x, 0) = —k and k = for t 7^ 0. We find 

m -t)\ - -^^m^M-^rt) 

n(x, t)U.,. - 2^A;-i -exp(-Art)) ' ^^'^"^ 

n(x,-t)|,.p. = 0. (5.86) 

Hence, according to equation ()4.12j) . the survival probability in mean field approximation 
is given by 

-2Tlg" X 2T/g" 



Pit) = 1 - exp {- 



.exp(Art) — 1/ exp(Art) — 1 

2/{\g"t) for r = 0, 

i2T/g") exp(-Art) for r > 0, t ^ 00 , (5.9) 

(2 |r| //)(1 + exp(-A |r|t)) for r < 0, t ^ 00 . 

We see that P{t) owns all the required properties displayed in equations ()4.14[ l^!T^ . 
According to equation ()4.11|) the inverse Laplace transform of 

1 - p(k, t) = (exp(n(0, -t) - kj<f{0))) = exp n(0, -t) Up. (5.10) 

leads back to 

PNit) = {5{N -J\f{t))expn{0,0)) = — / dA; (1 - P(A;, t))e'^^ 

J (T—ioo 



1 



(T+ioo 



kN 



= Po{t)5{N) + ^ / dA; (P(cx), t) - P{k, t))e 

=: Po{t)6{N) + P^it) . (5.11) 

Po{t) = 1 — P(t) is the extinction probability, and we have used equation ()4.12|) . The 
subtraction in the integral that defines -P/v(t) is needed for the application of Jordan's 
lemma. Only after this subtraction the inverse Laplace transform can be calculated 
by closing the integration path with a large semi-circle which tends to infinity in the 
left complex /c-plane. By this subtraction-procedure the (5-contribution of the absorbing 
state is splitted off. Using the saddle point solution ()5.8ajl we find 

Po{t) =exp{-A{t)N{t)), (5.12a) 
PN{t) = A{t) sj N{t)/N h (2A{t) N{t)N^ exp[-A(t)(iV(t) + A^)] . (5.12&) 
Here we have defined 

m=-, ^TTT' (^-^^fl) 

1 — exp(— Arrj 

iV(t) = exp(-Art) . (5.13&) 

Note that with the identifications A = 1 and g" = 2 the mean field solution ()5.12ajl turns 
into the distribution calculated with the Fokker-Planck equation ()2.4aj) with xq = 1 as 
we have anticipated in section II. 
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Asymptotically, if A{t)N{t) ^ 1 and A{t)N is finite, we can expand equations 
()5.12aj) and obtain 

PN{t) ^ (1 - A{t)N{t))S{N) + A{tyN{t) exp{-A{t)N) . (5.14) 

This expression can also be derived if expn(O,0) is approximated by 1 + n(0,0) in 
equation ()5.1H) . 

For logarithmic corrections to the mean field results in the upper critical dimension 
rfc = 4 see EH- 



6. Multiplicative noise and comparison with other work 

Now we equip the Langevin equation ()H.1|1 with multiplicative noise 



C(x, t)C(x', V) = Xg"n{^, tfdi^ - x')5(t - t') (6.1) 

instead of the absorptive one, equation ()3.2|) . It is well known that in this case the 
Langevin equation (j3.1|) cannot represent a physical reaction-diffusion system fl^ I14j. 
At least diffusional noise must be included and the response functional changes to 

/d III 
dtd'^x |An(^A~^— + (r - V^) + - ^—nr^n - An(Vn)| . (6.2) 

The upper critical dimension dc decreases from 4 to 2. The naive scaling dimensions of 
the fields become d for n and for h. Hence, one expects that also non-Gaussian noise 
contributions becomes relevant. Indeed the field theory based on J\ equation ()6.2p . is 
not renormalizable as it stands. However, a quasicanonical transformation to "bosonic" 
fields 

a"*" = 1 + a = exp h , (6.3a) 
a =nexp(— n), (6.35) 

leads to the elimination of the diffusional noise term. The resulting action is known 
from the theory of branching and annihilating random walks The corresponding 

Langevin equation is formally equipped with imaginary multiplicative noise. 

It is instructive to apply the saddle point approximation to the response functional 
()6.2|) . This procedure is expected to be correct above the upper critical dimension. The 
appropriate saddle point equations are given by ()5.2aj) with the replacement of g" by 
g"n^ and an additional contribution arising from the diffusional noise. However, both 
terms are cancelled for homogeneous solutions with n = 0. We get the time-dependent 
solution 

n(x, — t)|s.p. = —k exp(— Art) . (6.4) 
which goes to infinity for k —>■ oo. Thus, for each finite time t we obtain from equation 

iH 

P(t) = 1 — lim expf;,(x, — t)|s.p. = 1. (6.5) 

fc— >oo 

Hence, the survival probability for processes equipped with multiplicative behaves 
trivially at least in the saddle point approximation. It is obvious that an expansion 
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of expn is not allowed. Indeed, it is appropriate to change to the bosonic variables 

Now we are in a position to compare our result for the survival probability with 
the corresponding expression given by MGT. As we have seen, the distribution function 
PN(t) ()5.11|1 consists of two parts 

PN{t) = PN{t) = ((5(iV-Ar(t))expn(0,0)) 

= {1 - P{t))6{N) + P{t)pN{t) (6.6) 

where the nonsingular distribution pNit) is normalized to dN piy(t) = 1. The 
survival probability P{t) is a monotonic decreasing function in the case of absorptive 
noise and equal to one in the case of multiplicative one. MGT have introduced the 
expression 

P(*^^^)(t, a) = {e{J\f{t) - a)h{0, 0)) (6.7) 

as a substitute of the survival probability. In the case of absorptive noise, where 
asymptotically the expansion of exp h is allowed, we get from equation ()6.6p 

{6{N-Af{t))n{0,0)) = P(t)[p^(t) - 6{N)] . (6.8) 

Hence, one sees immediately 

limP(^^^^)(t,a) = P(t) (6.9) 

in the large t and small r limit. The advantage of our expression ()4.1Hj) for the 
survival probability is its calculability in perturbation theory. This contrasts to the 
substitute ()6.7j) . in which the Heaviside-^-function is a insurmountable complication for 
any calculation. 

Things become worse for multiplicative noise. As we have shown, the expansion of 
expn is dubious, and the survival probability P(t) = 1 as long as t is finite. Even if the 
distribution pN(t) piles up near N = with pN{t) ~ A^'' for small A^ and a width B(t), 
the substitute behaves like 

P^^'^^\t, a) = l + 0{{a/B{t)Y+^) (6.10) 

and goes to zero exponentially if a ^ B(t). Hence, a is a dangerous variable, and it is 
not clear which value should be assigned to it. 

7. Other growth models with absorptive noise 

In this section we will briefly discuss the applicability of our expressions to other growth 
processes with absorbing states widely discussed in the literature. We present several 
variants of the original DP, for which the expressions ()4.7p4.1H|l for the probabilities are 
well applicable. 
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7.1. Tricritical directed percolation 

Tricritical directed percolation (TDP) can be understood as the limit of DP in which 
the renormalized coupling constant g' goes to zero. In this limit no growth limiting 
term shows up in the original mean field part of the equation of motion of DP ()3.1I3.3|) . 
We expect that formerly irrelevant couplings become now relevant. The nature of these 
terms depends on the physics of the growth process under consideration. The simplest 
case is the occurrence of a higher order effective coupling in the reaction rate. In TDP 
IHH] one assumes the form 

R[n] = \{t + ^-n + L'^ (7.1) 
that leads to an interaction part of the response functional ()3.6p 

VTDp[n, n] = ^ (^g'n + ^n' - g"fi) . (7.2) 

In the mean field approximation, the tricritical point which separates the continuous 
DP-transition from a discontinuous one is found at t = g' = 0. The noise term has the 
absorbing form (|3.2j) . However, the dual symmetry (|3.12j) is lost and we expect different 
exponents (3 ^ [3 = [3' . Indeed, a field-theoretic renormalization group calculation up 
to 2-loop order in an e-expansion below the upper critical dimension (ic = 3, e = 3 — (i, 
yields 

/? = 1-O.458e + O(e2), j3' = l + 0{e^). (7.3) 

The other exponents are 2 = 2 + 0.0086 e + 0{e'^) and v = l/2 + 0.0075 e + 0{e^). The 
crossover exponent = 1/2 — 0.012 e + 0{e^) describes the scaling of the crossover from 
TDP to DP. It can be easily seen that the mean-field calculations of the probabilities 
in section V lead to the same results for TDP as found for DP. 



7.2. Compact directed percolation 

The universality class of compact directed percolation (CDP) that includes e.g. the voter 
model jni EZI EHl EHl J describes another universality class that can be reached in the limit 
(7' — > of DP. In contrast to TDP, in CDP a limiting term in the deterministic part of 
the continuum Langevin equation of motion ()3.1l3.3p does not exist a priori. Thus, in 
the active phase with r < 0, the density grows to infinity corresponding to a state with 
a completely filled lattice. Renormalizing the density of the filled lattice to the value 
n = 1, we expect that formerly irrelevant couplings even in the noise become relevant. A 
special physical condition comes from the demand that the empty and the filled lattice 
are both absorbing states and the exchange n ^ 1 — r ^ — r is a symmetry. The 
patchy spatial structure of the CDP-states motivates an expansion of the equation of 
motion in the small quantity n{l — n). Hence, the relevant terms of the reaction rate of 
the Langevin equation ()3.1|3.3j) are given by 

R N(x,i) ^(X' = ^^(1 - ^(X' t)- (7.4) 
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The reaction rate is zero at the critical point r = 0. For r 7^ 0, the model has a bias 
which prefers the state n = or n = 1 depending on the sign of r. The noise in the 
empty or the filled lattice is zero. This enforces Gaussian fluctuations of the noise of 
the form 



C(x, t)C(x', f) = Xgn{^, - n(x, t))5(x - x')5(t - t'). (7.5) 
equations (j7.4l7.5p lead to the response functional (j3.6p of GDP 

JcDp[n,n]= y"dtrf'^x|An(A-^^-V^n + ™(l-n)-|nn(l-n))} (7.6) 

This response functional is invariant under the transformation n 1 — n, n ^ —n, 
T —T. For r = 0, J^cDpln^n] transforms under duality (j3.1!2|) to an action which 
describes annihilating random walks ^Ol E]. Note that the noise then becomes 
imaginary in the corresponding Langevin equation. Physically this results from the 
anticorrelations present in the annihilating walk process. Gonsequently, the notion of 
density variables are inappropriate in such a case. Instead the fields n, h represent 
annihilation and creation operators of a bosonic field theory. The scaling dimensions 
of n and n are not related, and are given trivially by the naive dimensions which are d 
and 0, respectively. The upper critical dimension is d^ = 2. The renormalization of the 
bias T is intimately related to the exactly computable renormalization of the coupling 
constant g. A short presentation of the field theory is given in the appendix. 
We note here the well known order parameter exponents of GDP for d < 2 

P = 0, P' = l. (7.7) 

The dynamical and the correlation length exponents are z = 2 and u = 1/d. The 
crossover exponent to DP is found to he (p = 2/d. 

For d > 2, the mean-field results for the probabilities derived for DP are applicable. 
The survival probability for clusters generated from a localized seed in an otherwise 
empty lattice at and below d = 2, as well as the percolation probability follows from 
our general formulas ()4.7|1 and ()4.13|) to 

P(t)=t-'^/2/(rt'^/'), PooOc|r|. (7.8) 
7.3. Conserved diffusing secondary field 

In analogy to the extension of Model A to Model G from the Halperin-Hohenberg 
catalog of fundamental field-theoretic models for critical dynamics, one can consider 
an extension of DP to DP-G In Model G a conserved diffusing density which is 
not critical itself is coupled to the relaxing dynamics of an order parameter. In the 
DP-G universality class, a conserved density c(x, t) is coupled to the growing density 
field n(x, t) of the DP-process. The relevant Langevin dynamics is given by 

dtn = XV^n- x(t + + fc^n + (7.9a) 
d^c = V\-fc + an) + C (7.9 &) 
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where the Langevin forces are Gaussian absorptive noise (j3.2j) and Gaussian diffusion 
noise 

C'(x,t)C'(x',f) = -7V2(5(x - x')(5(t - t'). (7.10) 

The first model of this type (with a = 0) was introduced by Kree et ah to deal with a 
catalytic poison which may damage a growing population. Ten years later, van Wijland 
et al. jlH] have shown that a special reaction-diffusion system {A + B ^ 2B, B ^ A) 
corresponds in general to DP-C with a cross-diffusion coefficient a 0. However, 
if 0"/ > 0, the renormalization flow leads finally to a violation of the stability bound 
2cr/ < g of the mean field part of equations ()7.9aj) . which is interpreted as the occurrence 
of a fluctuation induced first order transition. In this sense, the KSS model [a = 0) 
defines a tricritical point of the critical DP-C model with af < and with its own 
universality class. Note that the limit 7 — > leads to another universality class [IH] 
called the conserved lattice gas (CLG) or Manna class ^47]. The systems belonging 
to this universality class exhibit infinitely many absorbing states characterized by the 
different distributions of a conserved quantity. A recent renormalization group approach 

is dubious. It proposes an upper critical dimension = Q, whereas simulations 

130] clearly exhibit an upper critical dimension d^ = 4. 

As the result of special symmetries, the DP-C model has the exact exponents z = 2 
and u = 2/d below the upper critical dimension d^ = 4. The special KSS model 
is invariant under the duality transformation ()H.12j) extended by c(x, t) — > c(x, — t). 
Hence, for a = we have the equality 

P = P' = l-^ + Oie') (7.11) 
lo 

where the expansion results from a 1-loop calculation In general, for a 7^ 0, dual 
symmetry is lost. In this case it turns out, that n does not renormalize. Its anomalous 
dimension is zero whereas the anomalous dimension of n is nontrivial. One obtains |^ 
finally 

/3 = 1, /?' = 1 -0.1065 £ + 0(e^). (7.12) 

Also for the DP-C as well as the CLG class, the mean field calculation of the survival 
and percolation probability lead to the formulas known from DP and is expected to be 
correct above the upper critical dimension. 

7.4- Dynamic isotropic percolation 

Contrary to the growth models considered up to now, the dynamic isotropic percolation 
(dip) modelled by the general epidemic process (GEP) jSI] leads to memory terms in 
the equation of motion for the growing density. The debris (or the immunes in the 
language of the GEP) is generated by the agent, i.e. the active particles (or the infected 
individuals). It suppresses the growth process itself and leads finally to the extinction 
of the disease (n = in each finite spatial region). The spatial distribution of the debris 
in the stationary hmit is described by the usual static percolation statistics. Thus, it is 
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more appropriate to consider the density of the debris as an order parameter field. It is 
proportional to the time integral of the density of the active particles 



m(x,t) = A / dt'n{x,t'). (7.13) 

J — oo 

The relevant terms in the reaction rate of the Langevin equation ()3.H) are now given 

by 

^N(x,i) = A(r + (7'm(x,t)). (7.14) 
Therefore, the response functional is follows as j2ZI 

dtd'^x {ATi(A-i- + (r - V^) + g'm - y^)^}- (7.15) 

dip is not invariant under the transformation ()3.12|) . However, the variant of the duality 
transformation 

a~^h(x,t) ^ —am(x, —t) (7-16) 

with a = \/ g' / g" transforms the response functional ()7.15|) onto itself. This shows that 
the particle density of the dual dynamic percolation can be identified with the field 
of the debris of the primal dynamic percolation process. Relation ()7.16p leads to the 
identity [3 — vz = (3 between the exponents defined in ()3.13|) . and shows that one has to 
identify P = Pp, where /3p is the order parameter exponent defined in static percolation 
theory. It is easy to see that the mean field result ()5.9|) for the survival probability is 
applicable above six spatial dimensions also in the case of dynamic isotropic percolation. 
For logarithmic corrections in the upper critical dimension dc = Q see [S^ • Whereas the 
formulas ()4.7| I4.13|) are correct in general, one now has 

Poo = a' lim (m(0, oo)e^'^^(°)) = a^pst oc {rf . (7.17) 

instead of equation ()4.8p as a consequence of the symmetry ()7.16|) . Here pst is the 
stationary homogeneous density of the debris in the active percolating phase. 



7. 5. Tricritical dynamic isotropic percolation 

As a last example we consider briefiy the tricritical variant of dynamic isotropic 
percolation (TdIP) jHHj I54j . TdIP as a tricritical point followed by a first order transition 
can be founded on a generalization of the GEP. One introduces in the reaction scheme 
of the GEP weakened individuals as a further species. This species is weakened by an 
infected neighbor but not immune. Hence, a second infection in its neighborhood makes 
this weakened individuals more easy susceptible. In the consequence, this property 
introduces an instability that may lead to a discontinuous epidemic transition. The 
response functional ()7.15|) is extended to 

dtd'^x {An(A-i- + (r - V') + fm+ ^m' - (7.18) 
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The tricritical point is given in mean field theory by r = / = 0. The upper critical 

dimension is d = = 5. The symmetry ()7.16|) is lost. Hence, we find four independent 

critical exponents. In a e-expansion with e = 5 — d we obtained 

S 74 9 
/?=|-^^ + 0(e^), (3' = l--e + 0{e'). (7.19) 

The mean field result for the survival probability is applicable above five spatial 
dimensions also in this case. For logarithmic corrections in the upper critical dimension 
dc = 5 see [SlI 



8. Epilogue 

Survival and percolation probabilities are most important quantities in the theory and 
in the application of spreading. In the past, one has missed calculable expressions in 
the field theory of this phenomena. We have seen that it is possible to derive such 
expressions, and that their calculations are feasible in perturbation expansions. We 
have presented the results in a simple mean field approximation. Extension to higher 
order of perturbation series are possible without difficulties. We have presented results 

in [aiiinaini!. 

The outstanding role of the form of the absorptive noise is stressed. So called 
multiplicative noise would lead to trivial non-decaying survival probabilities. In the 
last chapter we have shortly considered some fundamental growth models equipped 
with absorptive noise that are representations of known universality classes of spreading 
phenomena. All considerations are applicable also to models with long-range spreading 

m 
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Appendix: Renormalized field theory of compact directed percolation 

It is well known that GDP offers many exactly calculable properties. The response 
functional of the GDP process with bias ()7.fj|l equipped with a further relevant term to 
describe the crossover to DP is 

J[n,n] = J dtd'^x ^Xn(^X'^^ - V^n + (r + r)n(l - n) + rn^ - ^nn{l - j | 

= J dtd'^x^^Xn(^X'^^-V'^n + rn + Tn{l-n)-^nn{l-n)^y (A.l) 

The new coupling term with r > suppresses now the density of the active particles 
in the formerly compact filled regions in contrast to the other terms. The deterministic 



Survival and percolation probabilities 



22 



term proportional to r and the noise proportional to g are effectively present only at the 
interface between empty and filled regions These interfaces are broken up if r > 
and the process crosses over to DP. Note that in a mean field picture the transition line 
between the absorbing and the active phase is given by Tc = —r. 

The renormalized field theory of the critical voter model below d = dc = 2, 
that is ()A.1|) with r = r = 0, is usually developed by mapping the model to the 
problem of annihilating random walks. The mapping consists of two steps. At first a 
duality transformation is applied, then the diagrams contributing to the relevant vertex 
functions jlHl SI] are summed up exactly. This can be done also directly for the field 
theory defined by the functional (jA.lj) . Including the parameter r in the propagator 
G{t, q, r) = 6{t) exp(— A(g^ + r)t), we find for the primitively diverging vertex functions 
^, where k and k denote the numbers of the amputated n- and n-legs of the irreducible 
diagrams, exactly 



12,2 = -212,1 = 2XgD, (A.2) 
ri^2 = -2ri,i + 2{iuj + Xq^ + Xr) = -2XtD. (A.3) 



Here, 



is the correction that arises from the summation of all contributing diagrams, the chains 
of "bubbles". The single bubble integral is given by 

/(.;, g, r) = A dt e"^-* J G{t, | + p, r)G{t, ^-p,r) 

T{l + e/2) ( g2 ^^.-,/2 



£(47r)'i/2 4 2A. 

where e = 2 — d. We renormalize g i>j g —>■ g = Zg, and choose the renormalization 
point RP by (r + g^/4 + iuj/2X) = yU^, where fi is the usual convenient inverse length 
scale. Furthermore, we define 

9D\rp = , , =:g = A,ufi' (A.6) 



1 + gI 



RP 



with = r(l + e/2)/ (Air)'^^'^. We obtain finally the renormalization constant 

Z = . (A.7) 

1 — u/e 



Equations ()A.2|) show that r renormalizes in the same way as g, whereas r needs no 
renormalization: 



T ^ T 



r. (A. 



With the help of this procedure all £-poles are eliminated from the perturbation series. 

In the same way one can study the correlation functions with the composite interface 
field If = n(l — n). An insertion of the field (p requires a further renormalization. It 
can be easily seen that this renormalization is given hj ip ^ ip = Z~^(p. Thus, all 
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renormalizations of the response functional ()A.1|) follow from the renormahzation of 
this interface field. The cumulant of k n-fields, k n-fields and m yj-fields, starting 
from uncorrelated initial conditions with n(x, t = 0) = uq is denoted by G^j^.j^. The 
renormahzation group equation (RGE) for the Greens functions reads 

{^^d^ + (3du + KTdr + mnj ^^.^^.^({x, t}, r, r, rio, m, A, yu) = 0, (A. 9) 

where the RG-functions result from the logarithmic derivative of the Z-factor, n = 
— 51nZ/(91nyu|p, holding bare parameters constant. We get 

/t = u, /3=(-£: + m)u. (A.IO) 

At the stable fixed point of the RGE, -u* = £ > 0, we arrive the scaling form 

G'.,^;^({x, t}, r, r, no, n,) = /^'^+™^G,_^^„({/x, iH}, tI-\ rl-\ n^, (A.ll) 

All the critical exponents follow from these equations. The DP critical line, the line 
between the absorbing and the active phase for r > 0, results as Tc oc —r'^l'^. From the 
asymptotic solution of the RGE directly in two spatial dimensions, one can infer 

G.,^;„^({x, 0, r, r, no, u) = f'XiirG,-,.^J{U, ft}, X{l)r\ rl'^ no, X{l)u), 

X{iy^ = l- ulnl. (A.12) 

This general scaling relation determines all the logarithmic corrections known for GDP 
in (i = 2. One gets, e.g., for the mean interface density at the critical point r = and 
with a non-universal time constant to, asymptotically 

((^(x,t)) oc no(l - no){lnt/to)-' + O((lnt/to)~'), (A.13) 

whereas (n(x, t)) = no = const, a relation well known from the exact solution of the 
two-dimensional lattice voter model [SSIEII- 
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